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We  consider  the  temperature  dependence  of  the  exciton  linewidth  in  single  InAs  self-assembled  quantum 
dots.  We  show  that  in  cases  where  etched  mesas  are  used  to  isolate  the  dots,  the  magnitude  of  the  linear 
temperature  coefhcient  and  its  dependence  on  mesa  size  are  described  well  by  exciton  scattering  by  acoustic 
phonons  whose  lifetimes  are  given  by  phonon  scattering  from  the  mesa  interfaces.  This  work  shows  that 
phonon  scattering  at  interfaces  and  surfaces  typically  present  in  quantum  dot  structures  can  be  important  in 
dynamical  processes  in  quantum  dot  systems. 
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The  interaction  of  excitons  with  phonons  in  semiconduc¬ 
tor  quantum  dots  determines  the  temperature  dependence  of 
the  optical  response  and  contributes  to  the  dephasing  of  the 
optically  excited  quantum  states  of  the  dots.  These  optical 
properties  are  of  current  interest  in  connection  with  proposed 
use  of  quantum  dots  in  implementations  of  quantum  comput¬ 
ing.  At  low  temperatures  the  photoluminescence  spectra  of 
single  dots  have  a  sharp  central  feature  and  a  broad 
background. 1-5  The  problem  of  phonon  effects  in  the  optical 
spectra  associated  with  bound  electronic  states  in  a  crystal 
was  first  treated  in  the  context  of  spectra  of  impurities6-9  and 
bound  excitons  in  bulk  semiconductors.9  In  terms  of  a  per¬ 
turbation  treatment  of  the  electron-phonon  interaction  in  the 
case  of  acoustic  phonons  the  central  peak  in  the  absorption 
or  emission  spectrum  is  associated  with  the  term  involving 
no  phonons  and  commonly  referred  to  as  the  zero-phonon 
line  (ZPL).  It  is  an  optical  analog  of  the  Mossbauer  line  in 
the  impurity  nuclei  transitions.8 

The  experimentally  observed  peak  has  a  width  that  was 
found  to  increase  linearly  with  temperature  at  T<  50  K  in 
many  single  dot  studies.1-4  The  inverse  of  the  ZPL  width  is 
related  to  the  rate  of  the  polarization  decay,  i.e.,  the  dephas¬ 
ing  time,  and  has  been  studied  for  single  dots  and  for  arrays 
of  dots  in  four-wave  mixing  experiments.2,10  In  recent  four- 
wave  mixing  experiments  on  arrays  of  dots  the  temperature 
dependence  of  the  linewidth  was  found  to  be  of  an  exponen¬ 
tial  activated  form  rather  than  linear.11,12  Two  different 
mechanisms  involving  off-diagonal  transitions  between 
quantum  dot  electronic  states  have  been  proposed  as  expla¬ 
nations  of  ZPL  broadening,  real  transitions  from  ground  to 
excited  states  in  larger  dots13  such  as  those  resulting  from 
interface  fluctuations  during  molecular  beam  epitaxy  (MBE) 
growth  and  virtual  transitions  giving  higher  order  terms  in 
exciton-phonon  interaction.14  Using  realistic  values  of  inter¬ 
action  parameters,  we  find  that  for  smaller  dots,  such  as  those 
resulting  from  strain  controlled  self-organized  growth  at  in¬ 
terfaces,  neither  of  these  mechanisms  explains  the  order  of 
magnitude  of  the  temperature  dependence  of  the  optical  line- 
width  observed  in  single  dots. 

In  this  work  we  consider  the  case  of  single  self-assembled 
dots,  and  in  particular  dots  in  etched  mesas  introduced  to 
isolate  the  single  dots.1,5  In  the  more  recent  of  these  studies 
the  luminescence  ZPL  width  was  obtained  as  a  function  of 
temperature  for  various  values  of  the  linear  size  of  the  mesa. 


PACS  number(s):  78.67.Hc,  78.55.Cr 

The  temperature  dependence  was  found  to  be  fitted  well  by  a 
straight  line  r(7’)  =  r(0)  +  aT  with  slope  a  monotonically  de¬ 
creasing  with  increasing  mesa  size.  In  this  work  we  define 
the  linewidth  as  a  full  width  at  half  maximum.  The  observed 
increase  of  the  low  temperature  linewidth  T(0)  with  decreas¬ 
ing  mesa  size  could  have  a  contribution  from  fluctuating 
charge  occupation  at  the  lateral  mesa  interfaces.  However, 
because  carriers  occupy  deep  traps  at  the  side  walls  with 
large  activation  energy,  they  would  contribute  to  weak  expo¬ 
nentially  activated  variation  rather  than  a  linear  T  depen¬ 
dence  of  the  linewidth.  Therefore  the  observed  mesa  size 
dependence  of  a  does  not  originate  in  fluctuations  of  inter¬ 
face  charges.5 

We  use  the  independent  boson  model  in  order  to  evaluate 
the  optical  spectra  of  localized  electrons  coupled  to  acoustic 
phonons.8,9  In  the  subspace  of  the  electron-hole  pairs  the 
interaction  can  be  written  in  terms  of  phonon  and  confined 
exciton  operators  a  and  B.  To  the  linear  order  in  phonon 
operators 

#int=  2  -M„m(q,o-)(aq  fr+fl!q  J,  (1) 

n,m  q 

where  Mnm  is  the  interaction  matrix  element.  The  summa¬ 
tions  are  over  the  confined  states  of  the  interacting  electron- 
hole  pair  in  the  quantum  dot  and  the  wave  vector  and  polar¬ 
ization  of  bulk  phonons.  If  only  the  diagonal  part  of  /7inl, 
which  we  shall  call  HdmV  is  retained,  one  obtains  an  exactly 
solvable  model.  This  model  was  used  in  Ref.  9  to  evaluate 
the  spectral  function  of  a  bound  exciton  in  a  bulk  semicon¬ 
ductor.  We  shall  treat  the  off-diagonal  part  as  a  perturbation. 
Following  Ref.  8,  we  shall  also  take  into  account  the  finite 
lifetime  of  the  phonons  and  explicitly  consider  various  pho¬ 
non  lifetime  limiting  processes. 

The  absorption  spectrum  can  be  obtained  as  a  Fourier 
transformation  of  a  one-particle  time-dependent  exciton 
Green’s  function.  This  is  equivalent  to  using  a  two-particle 
interband  correlation  function  in  the  Kubo  formula  for  linear 
response  in  terms  of  electron  operators15  and  gives  time  de¬ 
pendence  of  light  induced  polarization.  Using  the  method  of 
linked  cluster  expansion,15  one  can  easily  generalize  the 
treatment  in  Ref.  9  to  include  finite  phonon  lifetimes  in  the 
derivation  of  the  electronic  correlation  functions.  Assuming 
that  phonons  are  in  thermal  equilibrium,  the  phonon  correla¬ 
tion  functions  are 
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(dq(t)dq{0))  =  exp {i(oqt  -  yqt)nq,  (2) 

(dq(0)d+(t))  =  exp(/&y  -  yqt)(nq  +1),  1 3=  0, 

where  ?!q=[exp(^,aiq/A:B7’)- 1]_1  and  the  phonon  lifetime  rq 
=  h/  yq.  We  then  obtain  the  following  result  for  the  ground 
state  exciton  Green’s  function. 

rco 

G(  a)  =  -  i  dtei(w~a,°+is),ef(t)(s  —  +  0) ,  (3) 

Jo 

where  ftco0  is  the  energy  of  the  exciton  state 


f(t)  =  (iA-T/2)t/h-f(t), 


r/2  =  E  |Moo(^’tr)|2rq 


75“A - 2~(2«q  +  0. 


q,cr  **'  “qcr 
,|2r*2„2 


fit) = 2 


lMoo(q^)l2(*2w5o-+  2/^rqwq(r) 


q,(T  (^2tv+  rip 


,A2 


(4) 

(5) 


X((2nq  +  1)(1  -e  V/fi)  +  {[1  -  cos(wqo.f)](2nq  +  1) 

-  /sin(«qo/)}e~VM),  (6) 


and  the  temperature  dependence  enters  through  2nq+l 
=  coth(ha>q/ kBT). 

In  the  evaluation  of  the  exciton-phonon  interaction  matrix 
elements  we  consider  a  cylindrically  shaped  In^Ga^As 
quantum  dot  (a  disk)  of  height  L  and  base  radius  R.  The 
Coulomb  correlation  in  the  electron-hole  pair  is  included 
through  variational  wave  functions  for  the  ground  and  ex¬ 
cited  pair  states.  We  construct  exciton  states  out  of  one- 
particle  confined  states,  which  are  eigenfunctions  of  z  pro¬ 
jection  of  orbital  angular  momentum,  with  z  chosen  along 
the  growth  direction.  The  one  particle  confinement  potentials 
obtained  from  the  band  gap  discontinuities  Uh  i=e,h  can  be 
written  as  Vi(ri,zl)=Ui9{ri-R)+Ui0{\zi\-L/2)+ SV^r^z,) 
where  the  nonseparable  part  SVi{r,z)  =  Ui6{\z\-LI2)6{ri 
-R)  can  be  treated  as  a  perturbation. 

For  the  disklike  case  of  2 R>L  we  neglect  the 
lateral  tunneling.  The  anisotropy  of  the  valence 
band  should  be  taken  into  account:  mxy  A  mz.  For 
the  ground  state  of  the  electron-hole  pair  we  take 
%(re,rh,ze,Zh)  =  R)Jo(rhx0l/  R)  exp 

(— <a:0|re— rh|), where  functions  f„(z)  refer  to  vertical  confine¬ 
ment,  Jm{r )  are  Bessel  functions.  For  typical  cases  where  the 
hole  in-plane  mass  mxy  is  substantially  larger  than  the 
electron  mass,  we  take  the  following  variational  form 
for  the  degenerate  first  excited  state  1F±1(re,rh,ze,zh) 
=  A^iz^/Kz^Joivoi  /  R)  Ji(rhxu  /  R)  exp  (-  ai|re-rh|) 
Xexp(+/9h),  where  <ph  is  the  in-plane  angular  coordinate  of 
the  hole.  In  these  two  expressions  aQ  and  ax  are  variational 
parameters,  xmk  are  positive  zeros  of  Bessel  function  Jm{x). 

In  the  case  of  an  In0.6Ga0.4As/GaAs  dot  the  band  gap 
discontinuity  A Eg  ~  0.74  eV  and  we  took  40%-60%  split  be¬ 
tween  valence  and  conduction  band  edges  in  evaluation  of 
confining  potentials  Ue  and  Uh.  The  electron  mass  in  the  dot 


me= 0.0404  in  units  of  free-electron  mass  m0,  and  0.0665  in 
the  barrier  material.  The  hole  mass  parameters  in  the  dot  are 
fn;, ,=0.34,  mhxr= 0.09.  The  value  of  the  heavy  hole  mass  for 
the  in-plane  dispersion  was  obtained  from  the  analytical  so¬ 
lution  for  the  infinite  barrier  case  with  the  appropriate  values 
of  Luttinger  parameters.16  For  a  disk  with  dimensions  L 
=  30  A  and  /?  =  1 00  A  we  obtain  a0R  =  0.15,  axR= 0.14.  The 
average  distance  between  the  electron  and  the  hole  in  the 
ground  state  is  (|rf,-r/l|2)1/2  =  75  A.  The  excited  state  and  the 
ground  state  are  separated  by  37.7  meV.  The  effect  of  the 
Coulomb  interaction  is  to  shift  the  ground  and  excited  state 
energies  by  similar  amounts,  -12.6  meV  and  -11  meV,  re¬ 
spectively.  When  the  parameters  are  such  that  1  /  a0  i  the 

Coulomb  correlation  of  the  electron  and  hole  is  small  com¬ 
pared  to  their  confinement  by  the  band  discontinuities.  In 
such  a  case  one  can  use  one -particle  wave  functions  in  the 
estimate  of  interaction  matrix  elements. 

At  lower  temperatures  we  consider  the  interaction  of  elec¬ 
trons  and  holes  with  acoustic  phonons  only.  We  included  the 
band  deformation  potential  interaction,17  the  piezoelectric  in¬ 
teraction,  and  the  electron-phonon  interaction  produced  by 
the  phonon  induced  deformation  of  the  confining  potential, 
also  known  as  the  “ripple  mechanism.”111  A  formal  expansion 
of  cxp[/(f)]  in  Eq.  (3)  produces  a  spectrum  in  the  form  of  the 
sum  of  terms  that  involve  n  phonon  transitions  with  n 
=  0,1,2,....  At  low  temperature  the  spectrum  has  a  form  of 
ZPL  plus  a  broad  background  from  many-phonon  transitions. 
The  overall  line  shape  can  be  deduced  by  considering  the 
short  and  long  time  behavior  of  the  integrand  in  Eq.  (3).  If 
we  assumed  an  infinite  phonon  lifetime,  i.e.,  set  y=0  in  Eq. 
(2),  then  T  =  0  and  at  short  times  t  Re/(?)  'x-t2,  while  at  large 
t  Re /(f)  oc  1  /f.  In  such  a  case  the  zero-phonon  line  is  not 
broadened  by  the  diagonal  part  of  the  exciton-phonon  inter¬ 
actions. 

Before  evaluating  the  effect  of  the  finite  phonon  lifetime 
in  Eq.  (3)  we  consider  the  effects  of  the  off  diagonal  terms  in 
Hmt  in  Eq.  (1).  The  real  transitions  to  excited  states  can  con¬ 
tribute  to  exciton  dephasing  in  larger  disklike  dots.13  The 
order  of  magnitude  of  the  corresponding  contribution  to  the 
exciton  linewidth,  r(0— *  1),  can  be  obtained  from  the  ground 
to  excited  confined  exciton  state  with  one -phonon  absorp¬ 
tion.  For  the  R=10  nm  dot  we  obtained  that  at  /  =  1 00  K 
r(0^1)~O(l(r4  /xeV).  Recently  the  effect  of  virtual  tran¬ 
sitions  in  Hint  was  treated  by  deriving  an  effective  quadratic 
coupling  in  exciton-phonon  interaction.14  Using  the  realistic 
values  for  deformation  potentials,  we  estimated  the  corre¬ 
sponding  contributions  from  Ref.  14  to  be  an  order  of  mag¬ 
nitude  smaller  then  the  observed  temperature  dependent  part 
of  ZPL  width  in  single  dot  experiments.5  We  conclude  then 
that  the  contribution  of  the  off  diagonal  terms  in  Eq.  (1)  to 
the  experimentally  observed  linewidth  for  smaller  dots  is 
negligible.  It  was  suggested  in  Ref.  8  that  direct  anharmonic 
terms  should  be  included  in  Eq.  (1).  It  was  later  shown, 
however,  that  the  inclusion  of  terms  quadratic  in  phonon  op¬ 
erators  in  the  diagonal  part  of  the  exciton-phonon  interaction 
does  not  lead  to  the  broadening  of  the  zero-phonon  line  in 
case  of  infinite  phonon  lifetime.19 

When  the  phonon  linewidth  y  is  included  in  the  treatment, 
the  resulting  width  1  of  the  zero-phonon  peak  in  the  absorp- 
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FIG.  1.  Phonon  linewidth  induced  by  the  boundary  scattering, 
shown  as  a  function  of  wave  vector  for  lateral  and  growth 
directions. 

tion  spectrum  of  the  dot  is  given  in  Eq.  (5).  We  evaluated 
several  sources  of  the  finite  lifetime  of  acoustic  phonons: 
anharmonicity  induced  three -phonon  processes,20  scattering 
of  phonons  by  impurities  and  alloy  fluctuations,  and  bound¬ 
ary  phonon  scattering.21  In  cases  where  etched  mesas  are 
present,  we  find  that  well  below  room  temperature,  the 
acoustic  phonon  lifetimes  are  limited  primarily  by  scattering 
at  the  side  surfaces  of  the  mesa  for  smaller  mesas  and  at  the 
top  surface  of  the  structure  for  large  mesas.  We  also  find  that 
in  the  evaluation  of  Eq.  (5)  with  yq  given  by  the  boundary 
phonon  scattering,  almost  all  of  the  resulting  value  of  zero- 
phonon  peak  linewidth  comes  from  the  band  deformation 
potentials  which  we  took  here  to  be  -9.3  eV  for  the  conduc¬ 
tion  band  and  3  eV  for  the  valence  band.17 

We  obtain  the  phonon  lifetime  from  the  steady  state  Bolt¬ 
zmann  equation  for  phonons.  This  is  similar  to  the  evaluation 
of  boundary  scattering  in  the  thermal  conductivity.22  Results 
similar  to  those  given  here  are  obtained  by  evaluating  the 
eigenmodes  of  an  acoustic  waveguide  with  rough  surfaces. 
The  approach  commonly  used  in  the  theory  of  wave  scatter¬ 
ing  from  rough  surfaces  is  that  these  surfaces  are  ergodic.23 
In  the  case  of  sound  scattering  this  is  justified  if  the  phonon 
spectrum  is  continuous.24  Then  the  eigenfrequencies  acquire 
an  imaginary  part,  even  when  the  scattering  is  completely 
elastic.25'27  Using  the  relaxation  time  approximation,  we  ob¬ 
tained  in  the  case  of  completely  diffusive  surfaces  1  /  t° 
~2v  i  /  A  for  scattering  from  the  mesa  sides  and  1  It 
~vz/Az  for  scattering  from  the  top  surface  of  the  GaAs  top 
layer.  Here  v  ,  and  vz  are  the  in  plane  and  z  components  of 
the  phonon  velocity,  A  is  the  side  length  of  the  mesa  which 
was  assumed  to  have  a  square  cross  section,  and  Az  is  the 
thickness  of  the  top  layer.  For  example,  a  completely  diffu¬ 
sive  scattering  of  phonons  by  the  faces  of  0.4  /rm  mesa  gives 
y°~  10  /xcV.  To  compare  this  to  the  three -phonon  processes, 
we  can  consider  the  most  efficient  of  the  latter,  the  Landau- 
Rumer  process  of  relaxation  of  a  transverse  acoustic  mode 
into  thermal  longitudinal  modes20  giving  fi  /  t„ 
=  6g2a>q(kBT)4(l+2vT/3vL)/(TTpvTuL4h2),  which  we  evaluate 
at  q  =  0.2  nm'1  and  take  Griineisen  constant  g= 2.  This  gives 
about  0.5  peV  at  50  K  and  becomes  comparable  to  the 
boundary  scattering  phonon  linewidth  only  above  100  K. 
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Inverse  mesa  size  1/  A  (1/gm) 

FIG.  2.  Coefficients  of  the  linear  temperature  dependence  of  the 
ZPL  width  (a,  solid  line)  and  of  Huang-Rhys  factor  (a,  broken 
line),  shown  as  functions  of  the  inverse  mesa  lateral  size. 

Thus  the  anharmonic  phonon-phonon  interactions  are  not  im¬ 
portant  in  the  evaluation  of  Eq.  (3)  below  100  K. 

In  the  case  of  partially  diffusive  surfaces  we  use  a  specu¬ 
larity  parameter  defined  by  the  boundary  condition  for  the 
nonequilibrium  part  8n(q,r)  of  the  phonon  distribution 
function  at  a  particular  boundary,  e.g.,  Sn(q:,  A.) 
=p{q)8n{-qz,  Az)  for  the  scattering  from  the  top  horizontal 
surface  at  z= Az.  As  was  shown  in  Ref.  21,  Jq=y^i\ 
-p{q)\l\_\  +p{qj\  and  one  can  use  a  Gaussian  approximation 
for  p(q):  p(q)  =  cxp(-4iT81q2)  where  the  asperity  parameter  S 
is  the  average  height  of  the  surface  bumps.21 

In  the  evaluation  of  Eq.  (5)  we  assume  the  following  lin¬ 
ear  dimensions:  R=10  nm,  L  =  6  nm  for  the  dot.  A, =6  nm, 
asperity  parameters  S=  12  nm  for  the  mesa  sides  and  8Z 
=  1  nm  for  the  top  (epitaxial)  surface.  The  phonon  linewidth 
is  shown  in  Fig.  1  as  a  function  of  wave  vector  q  =  (Q,qz)  in 
two  directions  in  the  momentum  space  for  A =0.4  pm. 

The  width  of  ZPL,  T,  as  a  function  of  temperature  T  and 
mesa  size  A  was  obtained  from  Eq.  (5).  We  find  that  in 
evaluation  of  electron-phonon  matrix  elements  in  Eq.  (5)  al¬ 
most  all  off  the  value  of  the  width  of  the  zero  phonon  line 
comes  from  the  deformation  potential  interaction,  while  the 
ripple  mechanism  contributes  about  1%  and  the  piezoelectric 
interaction  contributes  much  less  than  that.  The  height  of 


Time  (ps) 

FIG.  3.  Optically  induced  polarization  amplitude  for  a  single 
quantum  dot  in  0.4  pm  mesa,  normalized  to  one  for  three  values  of 
temperature  and  two  values  of  zero  temperature  linewidth,  T(0) 
=  1  peN  (solid  lines)  and  T(0)=100  peW  (broken  lines). 
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FIG.  4.  Optical  absorption  for  a  single  quantum  dot  in  0.4  gm 
mesa,  normalized  to  one  for  three  values  of  temperature  and  cen¬ 
tered  at  the  ZPL  energy,  with  T(0)=  1  /zeV. 


ZPL  is  determined  by  a  Debye-Waller  factor  e~s  with  Huang- 
Rhys  factor7-9  S  defined  as  the  real  part  of  the  time  indepen¬ 
dent  term  in  Eq.  (6) 


q,a 


|M00(q,o-)|2ft2w?o. 

iy  2  2  ,  2\2 


(2nq  +  1). 


(7) 


At  temperatures  above  5  K  we  find  that  both  the  linewidth 
and  Huang-Rhys  factor  are  given  by  nearly  linear  functions 
of  T,  T(T)-T(())-  dP  and  S(T) — S(0)  =  crT,  and  we  evaluated 
the  linear  coefficients  a  and  u  as  functions  of  the  mesa  size 
A.  The  linewidth  coefficient  a  varies  approximately  linearly 
with  the  inverse  of  the  mesa  size  while  the  Huang-Rhys  co¬ 
efficient  cr  varies  only  slightly  with  the  mesa  size.  The  coef¬ 
ficients  are  shown  in  Fig.  2  as  functions  of  1  /  A. 

In  the  evaluation  of  the  total  linewidth  we  should  also 
include  the  sample  dependent  term  I '(())  which  has  been 
found  to  vary  from  1  /xeV  to  200  /xeV  for  single  dots  in 


mesas.1,5  The  time  dependence  of  optically  induced  polariza¬ 
tion  amplitude28  is  given  by  exp{  Rc[/(f)]}  and  is  shown  in 
Fig.  3  in  normalized  form  for  three  values  of  temperature, 
the  mesa  size  A=0.4  /xm,  and  for  two  values  of  T(0).  The 
optical  absorption  is  given  by  Im  G(co )  and  is  shown  in  Fig. 
4  normalized  to  one  for  the  same  values  of  T  and  T(0)  as  in 
Fig.  3.  The  small  polariton  shift  A  is  of  order  of  temperature 
dependent  part  of  the  linewidth.  The  emission  spectra  can  be 
obtained  easily  by  assuming  an  equilibrium  distribution  of 
photons  in  the  emitted  light.15  The  spectra  are  similar  to 
those  in  Fig.  4  for  lower  temperatures  except  for  the  appear¬ 
ance  of  a  small  asymmetry. 

Comparing  our  results  with  the  experiments  on  single  dots 
in  mesas  we  conclude  that  a  linear  temperature  dependence 
of  the  ZPL  results  mostly  from  deformation  potential  inter¬ 
action  of  confined  exciton  with  acoustic  phonons  that  are 
scattered  by  the  interfaces,  with  scattering  from  the  surfaces 
formed  in  lateral  fabrication  giving  the  dominant  contribu¬ 
tion  to  the  broadening  of  the  acoustic  phonons. 

We  also  point  out  that  even  for  those  quantum  dot  experi¬ 
ments  that  do  not  have  explicit  etched  mesas,  the  effects  of 
interface  induced  phonon  linewidths  discussed  here  contrib¬ 
ute  to  the  ZPL  broadening  of  the  exciton  due  to  phonon 
scattering  at  interfaces  and  surfaces  formed  during  growth 
and  fabrication.  In  such  cases  one  obtains  a  smaller  value  for 
the  coefficient  a  than  those  obtained  here,  which  is  consis¬ 
tent  with  the  corresponding  experiments.2  In  addition,  the 
scattering  of  phonons  at  interfaces,  which  are  typically 
present  in  quantum  dot  structures,  can  affect  a  wide  range  of 
phonon  induced  dynamical  processes  in  quantum  dot  sys¬ 
tems. 
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